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Abstract 

In the present paper we discuss a generation of large antimatter regions with sizes 
exceeding the critical surviving size. In the modern epoch domains with high anti¬ 
matter density evolve to single galaxies with a peculiar content of anti-helium and 
anti-deuterium. 


1 Introduction 

Astrophysical observations indicate that our local Universe is filled mostly with baryons. 
If the baryon charge was zero from the beginning the total baryon excess in modern 
epoch should be much smaller than the observed density [1]. It means that the baryon 
charge has been dominating in the early Universe. The charge could be distributed 
uniformly throughout the Universe or some antimatter domains could exist. 

Characteristic size and abundance of such antimatter islands depends on initial 
conditions and physical parameters of a model. Main observational signatures of the 
former values are gamma radiation and presence of antinuclei in cosmic rays. Indeed 
annihilations that take place at the borders between regions of matter and antimatter 
would contribute to a diffuse gamma-ray background. The distortion of CMB spectrum 
by released energy could be insignificant depends on the parameters of the antibaryon 
islands. 

A close contact of coexisting matter and antimatter in the early epochs is almost 
unavoidable [18]. Equal number of matter and antimatter islands signihcantly con¬ 
tribute to the diffuse gamma-ray background [16, 17]. The gamma-ray flux in 100 
MeV range caused by this kind of annihilation would be below the observable one only 
in the case when the characteristic size of domains exceeds 10^ Mpc. This fact requires 
baryon domination over the whole volume of the Universe. 
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Since antinuclei are very unlikely to be formed in the proton-antiproton (or proton- 
proton) collisions their presence in cosmic rays would be a strong indication of cosmic 
antimatter. Antinuclei could be naturally created inside the islands with high an- 
tibaryon density as a result of nucleosynthesis. Nowadays there are no evidences of 
antimatter domains existence [3, 2]. This does not exclude the case when the Universe 
is composed of matter with relatively small insertions of antimatter islands [17]. The 
latter may be revealed as antistars distributed within our Galaxy [19] or even distant 
galaxies. The absence of antihelium in the cosmic rays and annihilation signals leads 
to the conclusion that their fraction in the galaxy is smaller than 10“'^. Analysis per¬ 
formed in [22] indicates that antimatter islands must be separated from a space hlled 
with matter at least by the distance of about 1 Mpc. Nevertheless, some room even 
for antistars in our Galaxy still remains [24]. 

Formation of antimatter islands of different sizes in the Universe is the subject 
of many baryogenesis models [4]. They all are based on the assumption of explicit 
G- and GP- violation. A detailed review of hrst models can be found in [6], [7]. The 
models of spontaneous baryogenesis are described in [8], [9], [10], [11]. Another possible 
mechanism for antimatter domains origin is based on the class of Affleck-Dyne (AD) 
baryogenesis models. The foundation of AD baryogenesis is a peculiar evolution of 
the scalar field which carries the baryonic charge, on the background of cosmological 
expansion. This held is supposed to be coupled to supersymmetric particles, that store 
the baryon number when scalar held decays. The important feature of these models is 
the existence of ”hat directions” in held space, in which the scalar potential vanishes 
[10], [12], [13], [14], [15]. If some component of the complex scalar held moves along 
a hat direction, it can be considered as a free massless held, the so-called AD held. 
Usually this held behaves as Nambu-Goldstone (NG) boson. 

The main aim of our paper is to show that the antimatter domains can be manifested 
in the form of single galaxies with non-trivial chemical composition. 

The proposed approach was based on the mechanism of spontaneous baryogenesis, 
and necessarily implied a complex scalar held carrying the baryonic charge. Necessary 
attribute of the previous paper [20] is the presence of large amplitude of initial baryon 
number huctuations on the biggest cosmological scales, which leads to large amplitude 
of isocurvature huctuations on large scales what contradicts the GOBE data [21]. One 
of the purpose is to settle this issue and further develop the model. 

The paper is organized as follows. In Sec. II we establish the model and discuss 
the quantum behavior of the scalar held in the inhationary period. Section III contains 
the size distribution of antimatter domains. In Section IV we discuss the observational 
manifestations. In Section V we conclude. 


2 Setup 

In this section we present a main tools of the model [20] which is the basis of our 
study. Necessary corrections are also inserted to avoid the problem of the large scale 
huctuations. The essence of the problem is as follows. According to the model [20] 
antimatter domains of large scale are effectively produced due to fluctuations of some 
hypothetical scalar held. Unfortunately the amplitude of baryon fluctuations is too 
large what contradicts the GOBE data (and the Planck data as well). To get rid of 
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this problem one should suppress large amplitude fluctuations on the largest cosmo¬ 
logical scales. In this paper we show that the effective suppression can be achieved by 
modifying the original potential used in paper [20]. 

= +A"(l-cos(0)), (1) 

where 'I' = x(^)exp(i0), x(t) is a radial component of 'I' field with minimum at f jV^ 
and 0 is a massless NG field. The first term looks similar to the Higgs potential which 
is known to be modified by quantum corrections [23]. As the result it alternates its 
form depending on the top quark mass. Moreover landscape ideas [25], [15], [26], [27], 
[28] imply potentials of arbitrary forms. Here we use this freedom to propose specific 
form of the potential 


V{x) = Vo{x)-F{x), (2) 

suitable for our purposes. The modified potential (2) is close to the original one while 
X is small. 

During inflation, when the friction term is large the classical motion of the angular 
field 6 is frosen due to the smallness of the potential tilt. The dynamics of the radial 
field X is governed by the equation of motion 


df^ 


+ 3H 


dxW 

dt 


, dH(x) 
dx 


= 0 


( 3 ) 


where H is the Hubble parameter. 

While the field x is moving classically to the potential minimum x — flV^, the 
phase 6 = ipj f changes due to the quantum fluctuations against the de Sitter back¬ 
ground [30], [31]. During inflation the amplitude of the phase 9 changes as 


69 ~ 


H 

27 r • |x(t)l 


( 4 ) 


every e-fold, where |x(t)| is the average modulus of the field x during that e-fold. 

The amplitude of baryon fluctuations 6Bi is proportional to 69 of the phase 9 [36], 
[20]. In the mentioned paper [20] x(t) = / = const so that the amplitude fluctuations 
69 are of the same order on different scales. This restriction is crucial for fitting the 
observational data. In our case the radial field x(f) varies with time according to the 
classical equation (3) so that 69 could be small on large scales comparing to those on 
small scales. This is the main difference between our model and the basic one [20]. 

Fig. I shows evolution of 69 in the modified potential. As one can see the phase 
fluctuations 69 remain small from the beginning and start growing after the 50 e-fold 
when the largest scales have been formed. Isocurvature fluctuations contribute to the 
CMB anisotropy as 


6T IHb 


( 5 ) 
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where (5 b. is the amplitu(ie of the initial baryon number fluctuations and Oo(^b) are 
the total (baryon) density in units of critical density. To avoid contradictions with 
COBE data one needs to have 8Bi ^ 69 < 10“^ in first 6 — 8 e-folds after the beginning 
of the inflation. With the mechanism proposed we have 66 ~ 10“®. Hence the conflict 
with observations is settled. 

The Planck data are satisfied as well. Recent results [29] put a strong constraint 
on the isocurvature mode of fluctuations. The latter is characterized by the ratio [21] 


A, 


1 / ^b \ 

ISvr \ Ho / 



( 6 ) 


Since fluctuation arising in our mechanism is similar to the axion mode, we use 
constraints on axion isocurvature mode f3iso < 0.039 while Piso < 0.075 in general case 
for the GDI component (CDM and baryon isocurvature mode) at low wavenumbers. 
The value of tensor spectral index —riT ^ 1 and varies depending on an inflationary 
model. It is possible to satisfy constraints on /3iso by choosing —tit in range 10“^-^10~^. 

The baryon charge of the field ^ should be accumulated in matter helds. To make 
it possible, an interaction of the matter and the held 'k is postulated in the form 



Figure 1: Computed values of 69 as a function of (60-N), where is a number of e-folds, 
N = H ■ t; The parameter values are: / = SOFT, A = 8 ■ 10“®, C = 4 • a = 1.455, 

Xo = 8.6 ■ lO^H 


Lint = g'^QL -k h.c. (7) 

where helds Q and L represent a heavy quark and lepton helds coupled to the ordinary 
matter helds. We shall assume that helds Q and ^ possess baryon number, while held 
L does not. The U(l) symmetry that corresponds to the baryon number is expressed 
by following transformations: 
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'I' —> 'I'exp(i/3), 
Q Qex.p{i/3), 
L 


( 8 ) 

(9) 


After the end of inflation the 6 field starts to oscillate around zero with the frequency 


uj = m 


-1 

B 


A2 


( 10 ) 


Once oscillations begin, the energy density of the PNG field converts to baryons and 
antibaryons. The expression for baryon number 




Svr^ 





dw 



( 11 ) 


can be found in [9], [20], where fig. is a volume containing the phase value 9i and 
coupling constant <7 1. 


3 Size distribution of antimatter domains 


In this Section we discuss the distribution of the antimatter domains following [20]. 
The total volume of all antimatter domains formed at the e-fold before the end of 
inflation may be calculated using the recursive procedure: suppose the total volume of 
all the domains with the average phase 9 formed by that time is V{9, Nt). Then their 
total volume at the {Nt — 1)**^ e-fold before the end of the inflation is given by [20] 


VilNt 


l) = e^V{9,Nt) + 


Vu{Nt) 


e^V{9,Nt) 


P{9,Nt-l)-h. 


( 12 ) 


where Vu{Nt) is the volume of the universe Vu{Nt) ~ g3-(Nu-Nt) e-folds 
before the end of inflation; Nu ~ 60 is the total number of e-folds during inflation. 
P{9,Nt — 1) gives the Gaussian distribution of the phase [30]: 


P{0,Nt) 


-exp 


H 


(^Nt = 


27r-\x{Nt)\ 


VNu - Nt 


(13) 

(14) 


Note that here Nt is the number of e-folds before the end of inflation, so it decreases 
with time; so if the moment Nq = Njj ~ 60, (t = 0) corresponds to the beginning of 
inflation, while in the end of inflation Nr = 0. Accordingly, the volume of antimatter 
domains in the beginning of inflation is V{9, Njj) = 0. 

The first term in the equation (12) - e^V{9, Nt) - is the total volume of antimatter 
domains formed before the Nfl^ e-fold. The second term 


v{9,Nt) 


Vu{Nt)-e^V{9,Nt) 


P{9,Nt 


1) • h. 


(15) 
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is the total volume of antimatter domains formed during the Nt^'^ e-fold. As the initial 
volume of each domain is H~^, the number of domains formed during a given e-fold is 


n = 


v{9,Nt) 


(16) 


Domains grow in size during inflation, so the earliest domains to form become the 
biggest at the end of inflation. Their linear sizes at present time are defined by the 
equation 

L{Nt) = 6 • (17) 

Here L{Nt) is the size in Mpc of the antimatter domain formed at the Nt^'^ e-fold. 

The result of numerical computations of the antimatter domains spectrum is repre¬ 
sented in Fig. 2. As one can see the maximum size of the largescale antimatter domains 
is about 300 kpc in the present epoch. 



Figure 2 : Spectrum of antimatter domains, with parameters of the potential ( 2 ): / = 80i7, 
A = 8 • 10-^ (F = 4 ■ a = 1.455, Xo = 8.6 ■ 10^77, 60 = ^ 


4 Observational manifestations 

During the inflationary process antimatter domains of different size arise. If a domain 
survives till the temperature 100 keV the nucleosynthesis starts inside it. Antinuclei 
formation begins with the reaction 


p -|- n —^ d -|- 7 


(18) 
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This reaction strongly depends on an antibaryon density inside the domain. For the 
reaction to be effective following condition has to be satisfied 

> H, (19) 

where (crv) = 4.55 • 10“^° cm^/sec [32] and one can use a relation H = T'^/M*i to 
go from the Hubble parameter H to the temperature T. Hence the anti-deuterium 
production takes place if the antibaryon density satisfy the inequality 

r^nucl > ^q18 ^20) 

The value of depends on the parameters such as g, A, / and the value of 6 in 
the certain domain. Since this mechanism is responsible for baryon charge generation 
it has to produce appropriate amount of baryon asymmetry. The observable number 
density of baryons is characterized by the ratio 

Ab = = 0.86 • 10“^° (21) 

s 

where s - is the entropy density and An® - is the baryon excess. 

According to [20]: 



where Y{9) = fg ,Mpi = 10^H, 5 * = 106.75. It is easy to obtain the value 
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Ab = 0.86 • 10“^*^ varying parameters are g and A with selected values of / and Oq. 

There is a wide room in the g — A parameter space to create the conditions for 
nucleosynthesis to begin. Fig. 3 represents a g — A parameter space for the largest 
domains that appeared at 50*^ e-fold with the biggest primordial antibaryon density 
(since we are interested in the largescale antimatter domains that could possibly be 
conserved till present epoch). In the modern epoch the size of these domains is about 
300 kpc. 

The right area and the area in the bottom left corner Fig. 3 are inscribed as ”no 
nucleosynthesis” which means that with these values of g and A there will be no antideu¬ 
terium production (18) in the largest domains. Preferred values (with nucleosynthesis 
going and Ab = 0.86-10“^*^) of g and A are represented by the line. All the values that 
lie above the line will lead to the very small baryon asymmetry and the ones below the 
line lead to the excessive amount of baryon asymmetry. 

The same parameter space plots could be drawn for every antimatter domain. 
Apparently there are too many parameters to take into account all of them, but we 
can do slices of parameter space as one presented in the Fig. 3. Hereafter we will 
use g = 0.01 and A = 0.64 to have baryosynthesis in most of largescale domains and 
satisfy (21). Note that there exist domains of various sizes (up to 300 kpc) with 

10"^^ < r/B < 6 • 10"^°. 

Due to the phase fluctuation, antibaryon density inside some domains can be very 
high in comparison with the background density. High density of antibaryons leads 
to the nontrivial chemical content after nucleosynthesis. Table. 4 shows the number 
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Figure 3: Parameter space for the largest domains with biggest primordial antibaryon den¬ 
sity. 


of antimatter domains with densities exceeding the Universe average density. These 
regions contain more '^He and Li, but the fraction of D and ^He is reduced. 

Note that this picture is similar to the Inhomogeneous Big Bang Nucleosynthe¬ 
sis (IBBN) [39]. In our case baryons are replaced by antibaryons and the difference 
between densities in different regions is not huge (IBBN model allow tjb > !)■ 

Table. 4 shows the abundances of antinuclei formed in the high-density antimatter 
domains. We used AlterBBN code (Version 1.4) [41] for numerical calculations. 

Upper limit to the ratio rj^ is 1.1 • 10“® due to the BBN constraints. On the 
other hand no antinuclei are formed inside the regions with low antibaryon density and 
antiparticles remain ionized till recombination epoch. 

This mechanism can lead to various antibaryon densities inside the domains of the 
same size. The probability of low-density domain to appear is much higher than the 
one of high-density domain. Hence, low-density regions with ~ 0 are much more 
abundant so that high-density domains are mostly surrounded by a space with zero 
antibaryon charge. The size of these voids should not contradict CMB constraints. 

According to our numerical calculations the largest scale of low-density antimatter 
regions is about 0.3 Mpc at present and there is no conflict with CMB. That differs 
from so-called Patchwork Universe. Matter-antimatter borders experience the lack of 
particles what reduces the annihilation rate in the epoch of inevitable annihilation [17], 
[37]. 

After the recombination large-scale structures begin to form. According to [34], 








Table 1: Number of antimatter domains with very high densities. 


L, kpc \qB 

bo 

o 

1 

o 

9.4-10-^° 

1.1 ■ 10-® 

1.4-10-9 

1.6 - 10-9 

100 

1.4-10® 

4.2 ■ 10^ 

9.2 • 10® 

1.2 - 10® 

90 

37 

1.6 ■ 10^ 

193 

13 

0 

0 

14 

214 

6 

0 

0 

0 

5 

66 

0 

0 

0 

0 

2 

29 

0 

0 

0 

0 


Table 2: Abundances of elements 


Vb 


m/H 

^He/H 

^Li/H 

^Li/H 


7.8 - 10-®9 

2.497 - 10-1 

1.713 - 10-® 

8.894 - 10-9 

7.353 - 10-19 

7.625 - 10-1® 

7.158 - 10-19 

9.4-10-19 

2.514 - 10-1 

1.232 - 10-® 

8.050 - 10-9 

1.008 - 10-9 

5.564 - 10-1® 

9.935 - 10-19 

1.1 - 10-9 

2.528 - 10-1 

9.157- 10-9 

7.437 - 10-9 

1.273 - 10-9 

4.186 - 10-1® 

1.262 - 10-9 

1.4-10-9 

2.550 - 10-1 

5.542 - 10-9 

6.639 - 10-9 

1.727-10-9 

2.582 - 10-1® 

1.721 - 10-9 

1.6 - 10-9 

2.562 - 10-1 

4.065 - 10-9 

6.258 - 10-9 

1.999 - 10-9 

1.914 - 10-1® 

1.995 - 10-9 


high-density domains with masses exceeding the Jeans mass can evolve into compact 
stellar objects, globular clusters or even antigalaxies what also lower the annihilation 
signals. As a result we get compact antimatter object separated from the matter 
regions by low density voids. 

The condition when the domain decouple from the cosmological expansion is 


Mb > Mj 


4 ^ 

3 ^ y/p 


(23) 


where p is the background energy density and Vg is the velocity of sound. Antibaryon 
domains with masses 10^ 10® Mq are gravitationally unstable and could collapse 

to form bounded systems like GC or antigalaxies. Tabled represents the number of 
antibaryon domains and their mass. That kind of a globular cluster (maybe with 
peculiar chemical composition) can exist in the nearby regions of our Galaxy. 

Along with large-scale domains that remain practically unaffected by the annihi¬ 
lation process there also exists a great number of ’’small” domains that could not 
survive till present time. Any antimatter regions with a size less than critical survival 
size Lc ~ 1 kpc in the contemporary epoch must be eaten up by the annihilation 
process [32]. Nevertheless it is important to show that these regions do not imply a 
significant distortion in the GMB spectrum. 

Total amount of antibaryons within the domains must be small compared to the 
total baryon number of the Universe. As it follows from (11), number of antibaryons 
depends on initial value of phase 9 before the end of inflation. We suppose that the 
average step of phase is 56 = at each e-fold, so that the phase value is distributed 
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Table 3: Number of antibaryon domains and their mass 


N 

M/Mq 

511 

10^ 

2.3 • 10^ 

10« 

4.5 • 107 

10^ 

2 ■ 10« 

10^ 


non-uniformly. Regions filled with phase —vr < 9i < 0 acquire nonzero antibaryonic 
charge. At the end of the inflation a size of regions is about H~^e^, where N is a 
number of e-folds. Under the parameters we use, the total number of antibaryons in 
all annihilated regions is Ng ~ 10^® what is negligibly small in comparison with a total 
number of baryons in the Universe. After its annihilation there will be no significant 
contribution to CMB. 

Finally let us discuss the angular size of the antimatter domains and annihilation 
flux on Earth. With new WMAP9 data it is possible to resolve objects of angular size 
of about Aa ~ 0.1. To remain hidden large antimatter domains of order 30 300 

kpc should be located at distances such that their angular size is less than Aa ~ 0.1. 
The latter is estimated as Aa = L/l^ where L is a size of the domain, is an angular 
diameter distance. Thus to mask their presence antimatter domains should be placed 
at Id > 300 kpc depending on size {Id > 3 Mpc for the largest domains of the size of 
300 kpc ). 

Gamma flux near Earth is expressed as [32] 


_~ 2 6 • 10“^ I _ - _ 

dn ~ 10-7 cm-3 


dy -2 -1 -1 

— cm s sr 

IdJ 


(24) 


while the diffused photon background is given by [35] 


dd> 

dn 


~ 10 ® cm s ^ sr ^ 


(25) 


below 1 GeV at high latitudes. Therefore largest antimatter domain with highest 
antibaryon density should be positioned at the distance Id > 30 kpc to stay unseen. As 
the result the antimatter regions could be unseen even if they exist in form of globular 
clusters or diffused antimatter clouds near our Galaxy. 

Recent discovery of antiproton signal at the energies about 100 GeV made by both 
PAMELA [42] and AMS-02 [43] does not put a strong constraints on our model. The 
latter gives antiproton energies much smaller than 100 GeV and there is no realistic 
mechanism to increase it. Collisions of antiproton CR from nearby anti-galaxy with 
protons of galactic halos would result in multiple pion production due to reaction 
pp —)• 5 pions [44]. Energy loss of pions and their decay give rise to diffuse gamma 
radiation with E ~ 300 MeV. Such an energy is close to lower threshold of Eermi LAT 
energy range. 
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Another observational signature would be interaction of anti-micrometeors of extra- 
galactic origin with the Sun and the Moon. Anti-meteors of different mass can be 
formed within anti-stellar systems. Collision with the Sun or the Moon would lead to 
the signihcant release of energy 

E = 10^^ erg. (26) 

\lmgj 

Annihilation rate limits can be found in [45]. 


5 Conclusion 

In the present paper we have developed a mechanism for antimatter region formation 
elaborated in [20]. The shortcoming of the previous model is the prediction of too large 
fluctuations of baryonic density, which contradicts observational data. The mechanism 
we propose removes that contradiction. 

The model predicts the generation of large antimatter regions with sizes exceeding 
the critical surviving size. In the modern epoch domains with high antimatter density 
evolve to single galaxies or globular clusters with a high content of antihelium and 
antideuterium. Such galaxies are separated from ordinary matter by voids of size at 
least I Mpc, which is consistent with CMB data. 
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